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Preface 



Surfing on the Ocean of Numbers - why this title? Because this little book does not 
attempt to give theorems and rigorous proofs in the theory 7 of numbers. Instead it will 
attempt to throw r light on some properties of numbers, nota bene integers, through a 
study of the behaviour of large numbers of integers in order to draw 7 some reasonably 
certain conclusions or support already made conjectures. But no matter how far we 
extend our search or increase our samples in these studies we are in fact, in spite of 
more and more powerful technologies, merely skimming the surface of the immense 
sea of numbers. - Hence the title. 

Most books in Mathematics are used as reference books. I still consult my first 
Number Theory book which I bought in 1949 - Elementary Number Theory by 
Uspensky and Heaslet. It was when I was young and enthusiastic and dreamt about 
becoming a Mathematician. I am still enthusiastic but I became a Physicist instead. . 
However, I stayed on the theoretical side and avoided to have to much to do with 
things that can break. But even so experiments are a major source of knowledge and 
maybe this book shows a little of a Physicist’s approach to Mathematics. Most results 
are presented or supported by tables and graphs. All calculations have been carried out 
on a Pentium 100 Mhz laptop using Ubasic as a programming language. Finally, the 
author has tried to make a book which should be easy and pleasant to read. 

A word about the beauty 7 of Mathematics and Number Theory in particular. The 
crystallized truth of a theorem, where a whole spectrum of mathematical thoughts 
come together to form an entity; is like a painting where designs and colours merge 
mto a work of art. But sometimes it is not the finished result which is the most 
interesting - it could be the unsolved problem itself. Why? Maybe it is the challenge of 
getting somewhere with it or the hours and days of thinking and trying that occupy the 
mind in a positive sense different from the problems of our time. It all brings piece to 
the mind - it’s like walking in the silence of the forest enjoying the trees, the sun and 
the blue sky 7 , and should it happen that all the bits and pieces suddenly fall into place 
to give a solution then it is the most sublime experience for the human mind - eureka! 
But then the interest in the problem fades away unless solving the problem created 
new ones - and that is almost always the case. 
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Most topics in this book have been selected from Only Problems , Not Solutions by F. 
Smarandache. Others have bee suggested by Dr. R. Muller of Erhus University Press. 
A few problems which the author has found interesting originate from the Numbers 
Count Column of Personal Computer World. This journal has had great importance 
tor the author as a source of recreational Mathematics and I take this opportunity to 
thank the Editor of this column Mike Mudge for all correspondence and 
encouragement he gave me in the past. 

Illustrations, graphics, layout and final editing up to camera ready form has been done 
by the author. Tables have been created by direct transfer from computer tiles 
established at the time of computation to the manuscript so as to avoid typing errors. 

This book has come into being thanks to R. Muller at Erhus University Press who has 
never failed an opportunity to give his support and encouragement. Rapid e-mail 
exchange between him in the USA and me in France has greatly facilitated our work. I 
also thank Dr. Muller’s colleagues for their help. Many thanks are also due to my son 
Michael Ibstedt for his help and advice concerning computer equipment and software. 

Last but not least my warm thanks to my dear wife Anne-Marie for her encouragement 
and endless patience with a husband who does not always listen because his mind is 
somewhere else. 

February 1 997 
Henrv Ibstedt 
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Chapter I 



On Prime Numbers 



This chapter deals with some computational observations on prime numbers and their 
distribution. These computations only skim the surface of the ocean of integers but 
they give an idea of the general behavior of primes and often support some of the 
many conjectures that are made concerning primes. Computer programs have been 
written in Ubasic with extensive use of some of the built in functions of this language. 
In some cases the use of these functions will be illustrated with a few lmes of program 
code. Most results are given in tabular and/or graphical form. 

1. On the Sequence a-p„ + b 

In his book Only Problems , Not Solutions 1 F. 
question: 

If (a,b)=l, how r many primes does the progression 
ns {1,2,...}, contain* 7 

Already for a=l and b=2 we run mto the classical unsolved problem “Are there 
infinitely many twin primes ? ” The answer to how r mam - ? is certainly equally difficult 
for other sets of parameters a,b. However, some interesting information on how a-pn + 
b behaves will be obtained for the first 10,000.000 primes pn. 

Let m be the number of primes produced by a-pn + b for n<N, i.e. if a-pn + b is prime 
we can write a-pn + b = q ra where q m is prime. The following Ubasic program lmes 
have been used to determine whether a-pn + b is prime or not; 
while N< 10000001 
p=nxtprm(p) 
inc N 

c=a%*p+b% 

if nxtprm(c-l )=c then inc m 
wend 

The program has been implemented for a set of values of the parameters a and b. The 
result is shown in table 1. It is interesting to visualize the result. Because of the 
logarithmic behaviour of the distribution of primes it is reasonable represent m/N as a 



Smarandache asks the following 
a-pn + b, where pn is prime and 



1 Unsolved problem number 1 7. 
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function of logi 0 N rather than as a function of N. For this reason the value of m has 
been recorded during the computation for N=10 s iO 2 , 10 3 , 10 4 5 10 5 , 10 6 , and 10 7 . 



Table 1 . Number of primes m in the progression a-p r ^b for n<N 



a.b /N 


10 


102 


1G 3 


10* 


105 


106 


10 7 


1,2 


5 


25 


174 


1270 


10250 


36027 


738597 


2.1 


5 


25 


166 


1221 


9667 


82236 


711153 


3,2 


8 


47 


290 


2350 


18919 


160127 


1392733 


4,1 


3 


21 


145 


1108 


9314 


78676 


685069 


5.2 


5 


26 


188 


1492 


12020 


103010 


903165 


6,1 


7 


39 


277 


2175 


18019 


153925 


1342255 


7.2 


4 


23 


167 


1288 


10634 


91232 
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The number of primes m in the sequence a*pn+b for n = 1, 2, ... N is illustrated in 
figure 2 for a =1 and b = 2, 4, 6, 8, 10 and 12 where the ratio M=m/N is plotted 
against logN. The corresponding numerical results are given in table 2. 

The general appearance of the graphs for various values of the parameters (a,b) in 
apn+b is very similar. In particular figure lb shows an interesting picture of curves 
r unning parallel to one another and in particular to the one for (a,b)=(l ,2), that is the 
curve for p+2 which corresponds to prime twins for which we have the classical 
conjecture that there are infinitely many. This makes the following conclusion 
reasonable. 

Conjecture: The progression apn+b , (a,b)=l contains infinitely many prime 
numbers. 



2 . Prime Number Gaps 

Smarandache asked how r many primes there is in the progression apn+b. For a=l and 
b=2 the question is equivalent to ‘how many twin primes are there? 5 . Since we have a 
very stable conjecture that there are infinitely many we now' want to know something 
about their distribution and also about the distribution of other prime number gaps 
g-prw-pn- With a small change in the Ubasic program used in the previous section we 
can study the distribution of primes over gaps g=2, 4, 6, ... 

p=3 

while p<N 
q=p 

p=nxtprm(p) 

B-CMV2 

me ffu) ‘Count the number of gaps = p-q. 

if p(u)=0 then p(u)=q ‘Store the smallest prime for which the 

wend ‘ gap occurs in p(). 

This program was run for primes p<N=2 10 9 The result is shown in table 3, where f is 
the number of gaps g and p the prime number for which the gap first occurs, N<2- 10 7 . 
All gaps g<292 except 264, 278 and 286 are represented in the table which is 
arranged so that gaps gs2 (mod 6), g^4 (mod 6) and g=0 (mod 6) are found in 
separate columns . Gaps g=0 (mod 6) occur much more frequently than the other two. 
This is illustrated in figure 3 which also shows that In f as a function of g has a near 
linear behaviour for all three types. The “wild 55 behaviour fot gaps>250 would 
certainly correct itself if the range of primes in the study were extended. The area 
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below the curves for gs2 (mod 6) and g=4 (mod 6) are equal as will be shown shortly. 
The curve for g=2 (mod 6) behaves very well while the one for g=4 (mod 6) shows an 
interesting npple effect. In particular it shows a ''bump” for g=70 which showed up 
already m the smallest sample N<10° for which g=70 first appeared. What causes this 
high frequency for g=70? 

For a prime number p>5 we have ps±! (mod 6). Let q and p be two consecutive 
prunes forming the gap g=p-q. We distinguish between the following cases: 



Shift 

1. qsl (mod 6) and p^l (mod 6) => g^O (mod 6) ++ 

2. qsl (mod 6) and ps-1 (mod 6) => gs4 (mod 6) +- 

3 . qs-1 (mod 6) and p^I (mod 6) => g=2 (mod 6) -+ 

4. qs-1 (mod 6) andp^-1 (mod 6) => gsO (mod 6) - 



A sequence of consecutive primes (with the first prime =5) can be characterized by the 
shifts: 

-■}--+ ___ _ i i ! I — ■ ■+■+■ 4 4 *f t t ) t t ■ ■ ■ 



The longest sequence of consecutive primes = 1 (mod 6) for p< 10* is of length 18: 



450988159, 450988177, 450988207, 450988231, 450988241, 450988261, 
450988297, 450988333, 450988339, 450988381, 450988387, 450988399, 
450988411, 450988423, 450988441, 450988471, 450988477, 450988567 

and the longest sequence of consecutive primes = -1 (mod 6) for p<10 9 is of length 22: 

766319189, 766319201, 766319231, 766319237,766319249, 766319261, 
766319273. 766319291,766319339, 766319357,766319363, 766319369, 
766319423, 766319441, 766319453, 766319483,766319507, 766319549, 

766119573, 766319579, 766319621, 766319627 

Le; the numoer of 23ps=2, 4 or 0 (mod 6) be 12 , f 4 , and f respectively for p<N. Then 
we will have fh— fh if the last shift is t- otherwise we will have ii-fffi - 

In a simple model it would be reasonable to assume that, at an arbitrary point in the 
sequence of shifts, the probabilities of rinding the next shift to be ++, -+ or — are 
equal, i.e. if we define F^fs/ffc+friT-i). and F(pfb/(fo+f 2+£0 we would 

have F2=F4=0.25 and Fc=0.5. This is not the case. 
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Before looking into this let’s first consider a related question: Do primes=l (mod 6) 
(notation f+) occur with the same frequency as primes = -1 (mod 6) (notation f) ? 
Table 4 shows a study of the number of primes f. and f+ congruent to -1 respectively 1 
(mod 6) for primes less than 10* for k=l , 2, 3, ... 9. 

Within the range of this study we have f - > f +. However, the ratio r =(f . - fr)/(f - - f+) 
is decreasing. Will eventually f . < f + ? 

We have proved that f 2 = f 4 (assuming the last shift to be +-). We will now study the 
relative frequencies defined through 

F 2 =f 2 /(fo+f 2 +f4) F4=fV(fo+f 2 +fO F 0 =fo/(£>+f 2 +f 4 ) 

Again we have F 2 = F 4 and of course Fo=l - 2 F 2 . To study how F 2 varies as we 
increase the number of consecutive primes p<10 k the execution of the program for gap 
statistics was stopped for k= 1 , 2, ... 9 to produce the data shown in table 5. 

Table 4. Number of primes s -1 respectively 1 (mod 6). r=(f- - f*)/(f- + f+) 



k 


r 


2 


3 


4 


5 


6 


7 


8 


'.yfi 


f. 


i 


12 


85 


* 616 


4805 


39264 


332383 


2880936 


25424819 


U 


i 


11 


81 


611 


4785 


39232 


332194 


2880517 


25422713 


f--f+ 


0 


1 


4 


5 


20 


32 


189 


419 


2106 




2 


23 


166 


1227 


9590 


78496 


664577 


5761453 


50847532 


f*10* 


0 


435 


241 


41 


21 


4 


3 


0.7 


0.4 



Table 5. Prime number gap distribution (mod 6) for primes <10* 


k 


g22(mod6) 


gs4(mod6) 


gs0(mod6) 


Total 


F2 


1 


2 


1 


0 


3 


0.5 


2 


9 


8 


7 


24 


0.3541666 67 


3 


58 


56 


53 


167 


0.341317365 


4 


379 


378 


471 


1228 


0.308224756 


5 


2870 


2868 


3853 


9591 


0.299134605 


6 


22839 


22837 


32821 


78497 


0.290941055 


7 


189285 


189284 


286009 


664578 


0.284819088 


8 


1616471 


1616470 


2528513 


5761454 


0.280566416 


9 


14107250 


14107249 


22633034 


50847533 


0.277442162 
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Tobte3. 



f 


P 


9 s * 


f 


_P ... 


9=0 


f 





6388042 


3 


4 


6386967 


7 


6 


11407651 


23 


5069051 


89 


10 


6568071 


139 


12 


8472823 


199 


4690561 


113 


16 


3527160 


1831 


18 


6427670 


523 


3528810 


887 


22 


3030348 


1129 


24 


4600962 


1669 


2190452 


2477 


28 


2386944 


2971 


30 


4298663 


4297 


1359889 


5591 


34 


1430231 


1327 


36 


2341569 


9551 


1103677 


30593 


40 


1308406 


19333 


42 


1940894 


16141 


796213 


15683 


46 


687135 


81463 


48 


1190342 


28229 


678359 


31907 


52 


511183 


19609 


54 


856601 


35617 


444581 


82073 


58 


383239 


44293 


60 


789454 


43331 


247659 


34061 


64 


253846 


89689 


66 


466901 


162143 


191321 


134513 


70 


272834 


173359 


72 


277514 


31397 


137620 


404597 


76 


122523 


212701 


78 


233230 


188029 


119756 


542603 


82 


85030 


265621 


84 


176328 


461717 


63174 


155921 


88 


65612 


544279 


90 


133019 


404851 


44723 


927869 


94 


40821 


1100977 


96 


71864 


360653 


37946 


604073 


100 


39504 


396733 


102 


52752 


1444309 


24215 


1388483 


106 


20996 


1098847 


108 


36484 


2238823 


21894 


1468277 


112 


17316 


370261 


114 


26413 


492113 


11385 


5845193 


118 


10863 


1349533 


120 


23526 


1895359 


7408 


3117299 


124 


7521 


6752623 


126 


14443 


1671781 


5181 


3851459 


130 


7111 


5518687 


132 


8974 


1357201 


3881 


6958667 


136 


3380 


6371401 


138 


6567 


3826019 


3970 


7621259 


U2 


2393 


10343761 


144 


4104 


11981443 


1776 


6034247 


148 


1966 


2010733 


150 


4022 


136262S7 


1288 


8421251 


154 


1561 


4652353 


156 


2152 


17983717 


886 


49269581 


160 


1012 


33803689 


162 


1413 


39175217 


661 


20285099 


166 


553 


83751121 


168 


1271 


37305713 


607 


27915737 


172 


430 


38394127 


174 


729 


52721113 


332 


38089277 


178 


292 


39389989 


180 


638 


17051707 


238 


36271601 


184 


235 


79167733 


186 


342 


147684137 


124 


134065829 


190 


205 


142414669 


192 


219 


123454691 


109 


166726367 


196 


112 


70396393 


198 


221 


46006769 


91 


378043979 


202 


71 


107534587 


204 


129 


112098817 


44 


232423823 


208 


56 


192983851 


210 


141 


20831323 


35 


215949407 


214 


38 


253878403 


216 


50 


202551667 


21 


327966101 


220 


36 


47326693 


222 


31 


122164747 


18 


409866323 


226 


15 


519653371 


228 


21 


895858039 


17 


607010093 


232 


3 


525436489 


234 


23 


189695659 


10 


216668603 


238 


8 


673919143 


240 


15 


39199543! 


8 


367876529 


244 


5 


693103639 


246 


7 


555142061 


6 


191912783 


250 


8 


387096133 


252 


8 


630045137 


3 


1202442089 


256 


1 


1872851947 


258 


2 


1316355323 


3 


944192807 


262 


1 


1649328997 


270 


2 


1391048047 


1 


1438779821 


268 


1 


1579306789 


276 


1 


649580171 


1 


1851255191 


274 


1 


1282463269 


282 


3 


436273009 


2 


1667186459 


280 


2 


1855047163 


288 


2 


1294268491 


i 1 


1948819133 


- m. 


1 


1453168141 
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Rgure 3. In[f) as a function of g for N<2-10 9 . g^O (mod 6) upper solid line, gs2 (mod 
6) lower solid line and gM (mod 6) dashed line. 
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Figure 4. Relative frequency of prime gaps . F 2 = F* (thick line) and Fo (thin line). 

Conclusion: F# <0.45 and F 2 = F 4 >0.27 for primes <10*. 

The approach to the values 0.5 and 0.25 that one would have expected is very slow 
and is slowed down with increasing k, - one realizes that the interval 8<k<9 is ten 
times as large as the interval 0<k<8. The data used is cumulative but even if we 
consider only the interval between 10 8 and 10 9 we have F 2 = 12,659,767/45,680,669= 
0.2771 compared to F 2 =0.2774 for the whole interval between 0 and 10 9 . 

Question: Given an arbitrarily small number 8>0, does a prime pi exist so that 
Fr<0.25+5 for all p>pi? 
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Chapter II 



On Smarandache Functions 



1. Smarandache - Fibonacci Triplets 

We recall the definition of the Smarandache Function S(n): 

S(n) = the smallest positive integer such that S(n)! is divisible by n. 

and the Fibonacci recurrence formula: 



F b = F^i + F b -2 (n>2) 



which for F 0 = F i = 1 defines the Fibonacci senes. 



We will concern ourselves with isolated occurrences ot triplets n, n-1, n-2 for which 
S(n)=S(n-l)fS(n-2) and pose the questions: Are there infinitely many such triplets? Is 
there a method of finding such triplets which would indicate that there are in fact 
infinitely many of them? 

A straight forw ard search by applying the definition of the Smarandache Function to 
consecutive integers w us used to identify the first eleven triplets [1] which are listed 
m table 1 . As often in empirical number theory this merely scratched the surface of the 
ocean of integers. As can be seen from figure 1 the next triplet may occur for a value 
of n so large that a sequential search may be impractical and will not make us much 
wiser. 



Table 1 . The first 1 1 Smarandache-Rbonocci triplets 



# 


n 


Sin) 




Sfn-2! 


1 


11 


1 1 


5 


2-3 


2 


121 


2-11 


5 


17 


3 


4902 


43 


29 


2-7 


4 


26245 


181 


18 


163 


5 


32112 


223 


197 


2-13 


6 


64010 


173 


2-23 


127 


7 


368140 


233 


2-41 


151 


8 


415664 


313 


2-73 


167 


9 


2091206 


269 


2-101 


67 


10 


2519648 


1109 


2-101 


907 


11 


4573053 


569 


2-53 


463 
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However, an interesting observation can be made from the triplets already found. 
Apart from n=26245 the Smarandache-Fibonacci Triplets have in common that one 
member is two times a prime number while the other two members are prime 
numbers. This observation leads to a method to search for Smarandache Fibonacci 
triplets in which the following two theorems play a role: 

I. If n = ab with (a,b) = 1 and S(a) < S(b) then S(n) = S(b). 

II. If n = p a where p is a prime and a < p then S(p a ) = p. 




Figure 1 The values ofn for which the first 1 1 Smarandache-Fibonacci triplets occur. 



The search for Smarandache-Fibonacci triplets will be restricted to integers which 
meet the following requirements: 



n = xp a with a<p and S(x)<ap 


(i) 


n-1 = yq b with b<q and S(y)<bq 


(2) 


n-2 = zr c with c<r and S(z)<cr 


( 3 ) 
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Table 2a. Smarandache-Rbonacci triplets 



# 


N 


S(N) 




SfN-1 } 




SJN-2) 




I 


1 


4 


4 


* 


3 




2 


* 


0 


2 


11 


11 




5 




6 


* 


0 


3 


121 


22 


* 


5 




17 




0 


4 


4902 


43 




29 




14 


* 


-4 


5 


32112 


223 




197 




26 


* 


-1 


6 


64010 


173 




46 


* 


127 




-1 


7 


368140 


233 




82 


* 


151 




-1 


8 


415664 


313 




167 




146 


* 


-8 


9 


2091206 


269 




202 


* 


67 




-1 


TO 


2519648 


1109 




202 


* 


907 




0 


11 


4573053 


569 




106 


* 


463 




-3 


12 


7783364 


2591 




202 


* 


2389 




0 


13 


79269727 


2861 




2719 




142 


* 


10 


14 


136193976 


3433 




554 


* 


2879 




-1 


15 


321022289 


7589 




178 


* 


7411 




5 


16 


445810543 


1714 


* 


761 




953 




-1 


17 


559199345 


1129 




662 


* 


467 




-5 


18 


670994143 


6491 




838 


* 


5653 




-1 


19 


836250239 


9859 




482 


* 


9377 




1 


20 


893950202 


2213 




2062 


* 


151 




0 


21 


937203749 


10501 




10223 




278 


* 


-9 


22 


1041478032 


2647 




1286 


* 


1361 




-1 


23 


1148788154 


2467 




746 


* 


1721 




3 


24 


1305978672 


5653 




1514 


* 


4139 




0 


25 


1834527185 


3671 




634 


* 


3037 




-5 


26 


2390706171 


6661 




2642 


* 


4019 




0 


27 


2502250627 


2861 




2578 


* 


283 




-1 


28 


3969415464 


5801 




1198 


* 


4603 




-2 


29 


3970638169 


2066 


* 


643 




1423 




-6 


30 


4652535626 


3506 


* 


3307 




199 




0 


31 


6079276799 


3394 


* 


2837 




557 




-1 


32 


6493607750 


3049 




1262 


* 


1787 




5 


33 


6964546435 


2161 




1814 


* 


347 




-4 


34 


11329931930 


3023 




2026 


* 


997 




-4 


35 


11695098243 


12821 




1294 


* 


11527 




2 


36 


11777879792 


2174 


* 


1597 




577 




6 



21 






Table 2b. Smarandache-Rbonacci triplets 



# 


N 


S(N) 




S(N-1) 




S(N-2) 


t 


37 


13429326313 


4778 


* 


1597 




3181 


1 


38 


13849559620 


6883 




2474 


* 


4409 


1 


39 


14298230970 


2038 


* 


1847 




191 


7 


40 


14988125477 


3209 




2986 


* 


223 


2 


41 


1 7560225226 


4241 




3118 


* 


1123 


-2 


42 


18704681856 


3046 


* 


1823 




1223 


4 


43 


23283250475 


4562 


* 


463 




4099 


-10 


44 


25184038673 


5582 


* 


1951 




3631 


-2 


45 


29795026777 


11278 


* 


8819 




2459 


0 


46 


69481145903 


6X1 




3722 


* 


2579 


3 


47 


107456166733 


10562 


* 


6043 




4519 


-1 


48 


107722646054 


8222 


* 


6673 




1549 


-1 


49 


122311664350 


20626 


* 


10463 




10163 


0 


50 


126460024832 


6917 




2578 


* 


4X9 


11 


51 


155205225351 


8317 




4034 


* 


42X 


-5 


52 


196209376292 


7246 


* 


3257 




3989 


-5 


53 


210621762776 


6914 


* 


1567 




5347 


11 


54 


211939749997 


16774 


* 


11273 




5X1 


0 


55 


344645609138 


7226 


* 


2X3 




4423 


9 


56 


484400122414 


16811 




12658 


* 


4153 


-1 


57 


533671822944 


21089 




18118 


* 


2971 


0 


58 


620317662021 


21929 




20X2 


* 


1627 


0 


59 


703403257356 


13147 




10874 


* 


2273 


-2 


60 


859525157632 


14158 


* 


3557 




10601 


-5 


61 


898606860813 


19973 




13402 


* 


6571 


1 


62 


972733721905 


10267 




10214 


* 


53 


-4 


63 


1 185892343342 


18251 




12022 


* 


6229 


-2 


64 


1225392079121 


12202 


* 


9293 




2909 


-4 


65 


1294530625810 


17614 


* 


5X7 




11X7 


-3 


66 


1517767218627 


11617 




8318 


* 


3299 


-8 


67 


1905302845042 


22079 




21478 


* 


601 


-1 


68 


2679220490034 


11402 


* 


7459 




3943 


11 


69 


3043063820555 


14951 




12202 


* 


2749 


5 


70 


6098616817142 


24767 




20206 


* 


4561 


2 


71 


6505091986039 


31729 




19862 


* 


11867 


2 


72 


13666465868293 


28099 




16442 


* 


11657 


7 
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p,q and r are primes. We then have S(n)=ap,S(n-l)-bq and S(n-2)-cr. From this and 
by subtracting (2) from ( 1 ) and (3) from (2) we get 



ap = bq + cr 


(4) 


xp a - yq b = 1 


(5) 


vq b - zr c = 1 


(6) 



The greatest common divisor (p d ,q b )=l obviously divides the right hand side ot (5). 
This is the condition for (5) to have infinitely many solutions for each solution (p,q) to 
(4). Such solutions are found using Euclid's algorithm and can be written in the form: 

x = Xo + q b t, v = vo - p*t (5 ) 

where t is an integer and (xo, v 0 ) is the principal solution. 

Solutions to (5') are substituted m (6') in order to obtain integer solutions for z. 



z = (yq b - iyr c (6*) 

Solutions were generated for (a,b,c)=(2,l,l). (a,b s c)=( 1,2,1) and (a,b,c)=(l,l,2) with 
the parameter t restricted to the interval -ll<t<ll. The result is shown m table 2. 
Since the correctness of these calculations are easily verified from factorizations ot 
S(n), S(n-1) and S(n-2) these are given in table 3 for two large solutions taken from an 
extension of table 2. 



Table 3. Factorization of two Smarandache-Ebonacci triplets 



n = 
n-1 = 
n-2 = 


1 6,738,688,950,356=2 2 -3-31 • 1 93-15.269? 

1 6 738,688,950,355=5- 197-1 ,399- 1 ,74 1 -6,977 
1 6, 738,688,950,354=2- 7 2 - 19-23-53-31 3-23,561 


S(n)= 

S(n-1)= 

S(n-2)= 


2-15,269 
6,9 77 
23.561 


n = 
rvl = 
n-2 = 


19 448 047 080.036=2 2 -3 2 -43 2 - 1 7.093 2 
19 448 047 080 035=5-7-19-37-61-761-17,027 
1 9 448 047 080 034=2-97- 1 .609-3.631 • 1 7. 1 59 


S(n)= 

S(rvl)= 

S(n-2)= 


2-17,093 
1 7.027 
17.159 



Conjecture: 

There are infinitely many triplets n, n-1, n-2 such that S(n)=S(n-l)+S(n-2). 
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Questions: 



1 . It is interesting to note that there are only 7 cases in table 2 where S(n-2) is two 
times a prime number and that they all occur for relatively small values of n. 
Which is the next case? 

2. The solution for n=26245 stands out as a very interesting one Is it a unique case 
or is it a member of a family of Smarandache-Fibonacci triplets different from 
those studied here? 

References: 

C. Ashbacher and M. Mudge, Personal Computer World, October 1995, page 302. 



2. Radu’s Problem 



For a positive integer n, the Smarandache function S(n) is defined as the smallest 
positive integer such that S(n)! is divisible by n. Radu [1} noticed that for nearly all 
values of n up to 4800 there is always at least one prime number between S(n) and 
S(n+1) including possibly S(n) and S(n+1). The exceptions are n=224 for which 
S(n)=8 and S(n+l)=10 and n-2057 for which S(n)=22 and S(n+1)=2L Radu 
conjectured that, except for a finite set of numbers, there exists at least one prime 
number between S(n) and S(n+1). The conjecture does not hold if there are infinitely 
mam solutions to the following problem. 

Find consecutive integers n and n^l for which two consecutive primes pk and pk+i 
exist so that pk < Min(S(n),S(n~l )) and pt-i > Max(Sfn),S(n~- 1)). 



Consider 


D+l = Xpr S 


(1) 


and 


n = yprH 5 


(2) 



where pr and p^i are consecutive prime numbers. Subtract (2) from (1 ). 

XPr S - VPr-l S = 1 



(3) 
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The greatest common divisor (pr s ,|Vi s ) = 1 divides the right hand side of (3) which is 
die condition for this diophantme equation to have infmtely many solutions. We are 
interested in positive integer solutions (x,y) such that the following conditions are 
met. 



S(irH) = spr, i.e. S(x) < spt 


(4) 


S(n) = sjv-i , i.e. S(y) <spr-i 


(5) 



In addition we require that the interval 

spr s < q < spr-i s is prime free, i.e. that q is not a prime. 

Euclid’s algorithm is used to obtain principal solutions (xo, y 0 ) to (3). The general set 
of solutions to (3) is given by 

X=XO+Pr-l S , y = yo - pA (7) 



with t an integer. 

These algorithms were implemented for different values of the parameters d=pr-i-pr, s 
and t. The result was a very* large number of solutions. Table 4 show's the 20 smallest 
(in respect of n) solutions found There is no indication that the set would be finite. 
One pair of primes may produce several solutions. 

Wi thin the limits set by the design of the program the largest prime difference for 
which a solution was found was d=42 and the largest exponent which produced 
solutions was s=4. Some numerically large examples illustrating these facts are given 
in table 5. 



To see the relation between these large numbers and the corresponding values of the 
Smarandache function in table 5 the factorizations of these numbers are given below: 



11 822936647 1 52295784830 1 8 = 2-3-89* 1 93-43 ]• 1 6 1 278 1 2 

1 182293664715229578483017 = 509-3253-16128232 

1 1 157906497858100263738683634 = 2-7-372-56671-553333 

1 1 1 57906497858 100263738683635 = 3*5- 1 M9*- 1 6433-553373 

17549865213221 162413502236227 = 3-1 l 2 -307- 12671 -553333 

17549865213221162413502236226 = 2-23-37*71-419-743-553373 

2703299 7592 1 20525363470705 1 8228485 7039 1 3 1 4 = 2-3M7- 1 289-20 1 7- 1 1 9983- 167441* 

27032997592 1 20525363470705 1 82284857039 1313 = 37-231 17-2451 7-38303- 1 67443* 
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Table 4. The 20 smallest solutions which occurred for s=2 and d=2 



# 


n 


S(n) 


S(n+1) 


P* 


P2 


f 


1 


265225 


206 


202 


199 


211 


0 


2 


843637 


302 


298 


293 


307 


0 


3 


6530355 


122 


118 


113 


127 


-1 


4 


24652435 


926 


922 


919 


929 


0 


5 


35558770 


1046 


1042 


1039 


1049 


0 


6 


40201975 


142 


146 


139 


149 


1 


7 


45388758 


122 


118 


113 


127 


-4 


8 


46297822 


1142 


1138 


1129 


1151 


0 


9 


67697937 


214 


218 


211 


223 


0 


10 


138852445 


1646 


1642 


1637 


1657 


0 


11 


157906534 


1718 


1714 


1709 


1721 


0 


12 


171531580 


1766 


1762 


1759 


1777 


0 


13 


299441785 


2126 


2122 


2113 


2129 


0 


-4 


551787925 


2606 


2602 


2593 


2609 


0 


15 


1223918824 


3398 


3394 


3391 


3407 


0 


16 


1276553470 


3446 


3442 


3433 


3449 


0 


17 


1655870629 


3758 


3754 


3739 


3761 


0 


18 


1853717287 


3902 


3898 


3889 


3907 


0 


19 


1994004499 


3998 


3994 


3989 


4001 


0 


20 


2256222280 


4166 


4162 


4159 


4177 


0 


Table 5. Four numerically large solutions 



Pairs of consecutive integers 


■Ell 


d 


s 


t 




1 182293664715229578483018 
1 182293664715229578483017 


3225562 

3225646 


42 


2 


-2 


1612781 

1612823 


1 1 157906497858100263738683634 
1 1 157906497858100263738683635 


165999 

166011 


4 


3 


0 


55333 

55337 


17549865213221162413502236227 
1 7549865213221162413502236226 


165999 

166011 


4 


3 


-1 


55333 

55337 


270329975921205253634707051822848570391314 

270329975921205253634707051822848570391313 


669764 

669772 


2 


4 


0 


167441 

167443 



It is also interesting to see which are the nearest smaller P* and nearest bigger Pk-i 
primes to Si = Min(S(nXS(n+l)) and S 2 = Max(S(n),S(n+l)) respectively. This is 

shown in table 6 for the above examples. 



26 












